Quantum discord in a spin system with symmetry breaking 
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We analyze the quantum discord Q throughout the low-temperature phase diagram of the quantum 
XY model in transverse field. We first focus on the T — order-disorder quantum phase transition 
both in the symmetric ground state and in the symmetry broken one. Besides it, we highlight how 
Q displays clear anomalies also at a non critical value of the control parameter inside the ordered 
phase, where the ground state is completely factorized. We evidence how the phenomenon is in fact 
of collective nature and displays universal features. We also study Q at finite temperature. We show 
that, close to the quantum phase transition, Q exhibits quantum-classical crossover of the system 
with universal scaling behavior. We evidence a non trivial pattern of thermal correlations resulting 
from the factorization phenomenon. 



I. INTRODUCTION 

Correlations provide a characterization of many-body 
systems [T. In the quantum realm both classical and 
non local quantum correlations (like entanglement) are 
relevant. Although entanglement completely describes 
quantum correlations for pure states, it is in general sub- 
tler to characterize the pattern of correlations for mixed 
states. Indeed the quantitative interplay between classi- 
cal and quantum correlations has been formulated only 
recently with the introduction of the quantum discord, 
operatively defining quantum correlations in composite 
systems [SilSj. In this context, intense research activity 
has been devoted to spin models, generically displaying 
non trivial patterns of correlations in different interesting 
regimes [4 10 . 

Two peculiarities are relevant here: on one hand, the 
order-disorder Quantum Phase Transitions (QPT) occur- 
ring at zero temperature as long as the control parame- 
ter h is tuned across a critical value he [H]- It is worth 
noting that the quantum order arises because superse- 
lection rules lead to a symmetry breaking [T2]. Besides 
QPT, spin systems may exhibit a further remarkable phe- 
nomenon occurring at h = hf located within the ordered 
symmetry-broken phase, where the ground state is ex- 
actly factorized [CTll4j . and therefore correlations are ex- 
clusively classical. Such factorization consists in a quan- 
tum transition exclusively for entanglement [15] and is 
rigorously not accompanied by any change of symmetry. 

In this article we provide the details of the analysis of 
the quantum discord arising in the quantum XY spin- 
1/2 system, both at zero and at finite temperature [IB] . 
In particular we consider the ground state with broken 
symmetry. We analyze the quantum discord both at the 
QPT and close to hf. We evidence how the quantum crit- 
icality affects the quantum discord at low-temperature, 
allowing to define the quantum cross-over phenomena in 
an operative way. Close to the factorizing point, the 
finite-temperature quantum discord displays characteris- 
tic traits. 



The structure of this paper is as follows. In Sec.|TT]we 
give an overview on the notions of quantum and classical 
correlations in a general quantum system. In Sec. |III| we 
introduce a many body system that is suitable for our 
type of analysis: the XY model in an external trans- 
verse field. We also present some features related to its 
physics and to the way in which correlations are evalu- 
ated (Sec. IV). The analysis of our results is carried out 
in the subsequent sections, first for the zero-temperature 
case (Secjvl) and then when the temperature is switched 
on (Sec.|Vip. Finally in Sec. VII we draw our conclusions. 



II. 



QUANTUM, CLASSICAL AND TOTAL 
CORRELATIONS 



In a bipartite quantum system AB the total amount 
of correlations between A and B is given by the mutual 
information 



I{A:B)^Sip^)+S{p^)-S{p^^) 



(1) 



where S{p) — — Tr[/51og2 p] is the von Neumann entropy. 
In classical information, using the Bayes rule, an equiv- 
alent formulation of the mutual information is 



J{A:B) = S{A)-S{A\B), 



(2) 



where the conditional entropy S{A\B) = S{AB) — S{B) 
quantifies the ignorance on part A once a measurement 
on B is performed. In the quantum realm, a measurement 
generally perturbs the system and part of the informa- 
tion itself is lost. Therefore, when wc consider a quantum 
composite system, Eq. ([2| differs from Eq. ([T]). This dif- 
ference allows us to estimate the relative role of quantum 
and classical correlations in composite systems [5]. 

As a matter of fact, if we describe a measurement on 
part _B by a set of projectors {Bk}, then 



-AS 

P{k) 



— (i^0-Bfc)p'4S(i^^^^) (3) 
Pk 
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is the composite state conditioned to the fc-th outcome 
with probabihty pk = Tr[(iA «) Bk) p^^ (Ia ® Bk)]. This 
conditioned state is the key ingredient which allows to 
distinguish between classical and quantum correlations: 
in fact it generally differs from the pre-measurement state 
p^^, as well as the mutual information differs from classi- 
cal correlations. A reasonable definition of classical cor- 
relations C then comes by finding the set of measure- 
ments on {Bk} that disturbs the least the part A, i.e. by 
maximizing [21 [3] 

C{p^^) ^ max [S{p^) - S{p^^\{Bk})] . (4) 

{Bk} 

The difference between mutual information and classical 
correlations defines the so called quantum discord: 

Q[p^^)^I{A:B)-C{p^''). (5) 

In the estimate of quantum correlations between subsys- 
tems of a bipartite system, entanglement has been play- 
ing a leading role, in particular about the relevance of 
correlations in many body systems. However in general 
quantum discord differs from entanglement: for example, 
even if they are the same for pure states, they can display 
a very different behavior in mixed states. 

III. THE XY MODEL IN TRANSVERSE FIELD 

Hereafter we will focus on an interacting pair of spins 
1/2 in the anti- ferromagnetic XY chain with transverse 
field h. The Hamiltonian of the model 

describes the competition between two parts: the near- 
est neighbor anisotropic interaction on the xy plane, with 
an anisotropy tuned by varying 7 S (0, 1], and the cou- 
pling with external magnetic field h along the z direc- 
tion (throughout this paper we set h = ks — 1). Using 
a set of successive transformations ( Jordan- Wigner, Bo- 
goliubov, Fourier [H]), the Pauli matrices operators tr" 
(a = x, y, z) on sites j can be expressed in terms of oper- 
ators such that the Hamiltonian takes the diagonal form 

H = - ^ ^kvivk + const . (7) 

k 

The system is thus described as a gas of non interact- 
ing fermionic quasiparticles, where r/j, (77^) is the cre- 
ation (annihilation) operator of a fermion with momen- 
tum k. Furthermore, the Jordan- Wigner transformation 
allows an analytic expression for the correlation functions 
gaa{r) = (o'^ct"^.^) of any two spins in the chain far r sites 
from each other [18]. 

The exact solution of the XY model encouraged 
a plethora of studies concerning its critical phenom- 
ena [T51[Tn]. In particular, during the last decade, new 



insight has been made in the description of the physics of 
the system through the analysis of quantum correlations 
(i.e. entanglement) [Ij. 



A. The Phase Diagram 

The phase diagram of the XY model is marked by 
two peculiar values of the applied field h [H [TTl [13] . For 
7 G (0, 1] the system displays a zero-temperature (T = 0) 
continuum QPT at /ic = 1, of the Ising universality class 
with critical indices i^ = z = l,/3 = l/8 11J. In fact, 
for strong enough external fields [h ^ h^) all the spins 
tend to be aligned along the z direction, while the op- 
posite limit {h <C he) gives rise to a spontaneous mag- 
netization (with Z2 symmetry broken) along a direction 
on the xy plane, that is 7-dependent. At zero temper- 
ature, on the left side h < he oi the phase diagram the 
system is an ordered ferromagnet and the Z2 symmetry 
is broken, while on the right side h > he quantum fluctu- 
ations lead to a disordered phase where the system is a 
quantum paramagnet. At finite temperature T > the 
physics of the whole system is still affected by the QCP 
h — he at zero temperature. A V^-shaped diagram in the 
h — T plane emerges, characterized by the straight lines 
T = \h ~ he\ that mark the crossover region between the 
so called quantum critical region (T > \h — he\) and the 
quasi- classical regions surrounding it |llj . 

Besides the QCP, there is another value of the trans- 
verse field that characterizes the zero-temperature phase 
diagram of the XY model. Indeed it has been found 
that, at zero temperature and at a certain anisotropy 7, 
ioT h = hf = ^1 — 7^ the ground state is exactly factor- 
ized [13]: 

i*Ss>-ni'/'j)' (8) 

IV'J) being a normalized single-site pure state. Therefore 
it emerges that, even though the system is in a phase 
with very strong quantum correlations, there is a "criti- 
cal" set of values h f (7) where the state is completely clas- 
sical. This strange occurrence, regarded as a paradox in 
the first place [13], seems to be strongly connected with 
the reshuffling of correlations among the system. In fact, 
a deep analysis on the behavior of entanglement has re- 
markably shed new light on the relevant physics involved 
on hf 14J. In particular it has been shown that, tuning 
the external field from h<hf to h>hf, the entan- 
glement pattern swaps from parallel to anti-parallel l5]. 
Furthermore it has been observed that at zero temper- 
ature the bipartite entanglement has a logarithmically 
divergent range at hf, together with the fact that at fi- 
nite temperature there is a whole region fanning out from 
hf where no pairwise entanglement survives |20j . 
This strongly suggests that, along these critical values 
of field and temperature, the behavior of entanglement, 
and in general of correlations, plays a pivotal role in the 
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physics involved and hence in our understanding of it. 
In particular it seems that the interplay of correlations 
when the field is tuned across hf is the only accessible 
way, so far, to tackle the puzzling physics that leads to 
the factorized state of Eq. ([s]). Here we show that the 
quantum discord allows a fine structure of the phase dia- 
gram around he, and, most remarkably, it displays a non 
trivial scaling law at the factorization field hf. 



IV. CLASSICAL AND QUANTUM 
CORRELATIONS IN THE XY MODEL 

In order to compute Qr between any two spins A and B 
at distance r along the chain, the key ingredients are the 
single-site density matrices p^,p^ and the two-site den- 
sity matrix of the composite subsystem p"^^ [see Eq.([5|]. 
Due to translational invariance along the chain, single- 
site density matrices are the same for any spin and they 
are given by 



P = P = 



1 



9x 



9x 
1 - 9z 



(9) 



where ga = {^a) are the local expectation values of the 
magnetization along the three different axes. On the 
other hand, the expression of p"^^ may be cumbersome. 
In fact, the general two-site reduced density matrix for a 
Hamiltonian model with global phase flip symmetry has 
the following form ^21,: 



Pr 



/ A a a F 

a B C b 

a C B b 

V F b b D 



(10) 



in the basis {|00), |01), |10), |11)}, where |0) and |1) are 
eigenstates of (because of translational invariance, 
this density matrix depends only on the distance r be- 
tween the two spins: pr = p^^)- The various en- 



tries in Eq. ( 10 ) are related to the two-point correlators 
gapir) = i^f^j+r) local magnetizations, ac- 

cording to the following: 



A = 1 + gz +gzz, 
D = 1 - +9zz, 
'Q = l- 9zz 

C = 9xx + 9yy: 
F — 9xx 9yy 

express the parity coefhcients, while 

a = ffx + 9xz J 

b = ffi - 9xz 



(11) 



(12) 



explicit the contribution from the symmetry breaking. 

As long as the system is in the .Z2-symmetric phase, 
matrix elements in "low case" are null (a = b = 0). 



The symmetry breaking manifest itself in a,b ^ |21j . 
In the former case, the remaining non vanishing "up- 
per case" entries in Eq. (11) can be evaluated analyti- 



cally [TH], therefore we could use a fully analytical ap- 
proach to compute the quantum discord in the so called 
thermal ground state [6]. In this state the system ap- 
proaches the ground state by lowering the temperature 
towards the limit T = 0; for this reason the symmetry is 
conserved and the state is not in the "true" degenerate 
ground state. In the latter case, the Z2 symmetry is lost 
and beside the spontaneous magnetization g^, also the 
non trivially computable gxzij") correlation appears |22) . 
In that case we resort to the numerical Density Matrix 
Renormalization Group (DMRG) method for finite sys- 
tems with open boundaries |23j. 

Once we have access to the density matrices through 
the correlation functions, we can compute the explicit 
form of the mutual information and the classical corre- 
lation in order to distill the amount of pure quantum 
correlations in Eq. ([5| . Since the reduced density matrix 
of the single spin is the same for any site, the mutual 
information is given by: 



I{A : B) = 2S{p^) log2 Xu , 



(13) 



iy=0 



where Aiy(r) are eigenvalues of p^^ ■ If one further spe- 
cializes to Z2-symmetric states, Eq. ( 13 1 can be read- 



ily evaluated by using the following two ingredients: 
i) the single-site density matrix in Ec^. (|9| turns out 
to be diagonal, therefore its von Neumann entropy is 
S{p^) = Sbinii^ + 9z)/'2], where Sunip) = -p\og2P - 
(1 — p) log2(l — p) is the binary Shannon entropy; ii) in 
terms of the correlation functions, the eigenvalues of pf^ 
turn out to be [61 : 



•^0/1 = ^ (1 + 5;^^ ± \/.9^ + i9xx - 9yyy) 



X2/3 — 



9z 



± 



9yy\ 



(14) 



Once the mutual information (13) is known in terms 



of the correlation functions, one needs to perform a suit- 
able measurement on spin B in order to compute clas- 
sical C{pr) and eventually quantum Q{pr) correlations, 
as stated in Eq. ([s]). Following a procedure similar to 
Refs. m [21], we use a set of projectors {Bk} as lo- 
cal measurements on the spin B. In particular, working 
on the computational basis {|0), |1)} in the Hilbert space 
71% associated to the spin B, our general set of projectors 
is 



{Bk = VIlkV^} , fc = G,l 



(15) 



where Ilk = \k){k\ is related to the basis vectors and 
V G U{2) gives the generalization to any type of projector 
on B. As suggested in Ref. _6_, V can be parametrized 
as follows: 



V = 



' cosf sinfe-*"^^ 
^sinfe*-^ -cosf ^ 



(16) 
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where 9 G [0,tt] and 4> G [0, 27r) are respectively the az- 
imuthal and polar axes of a qubit over the Bloch sphere 
After a measurement has been performed on B, 



0.71 0.715 



in nl 



the reduced density matrix P will be in one of the 

following states: 



Po/i 



(17) 



This expression gives the explicit dependence of the sys- 
tem A with respect to the projective measurement per- 
formed on the spin B. The coefficient qk,a — <Zfc.a(^, 0) in 
the expansion depends on the projectors used to perform 
the measure on B (see Ref. for the explicit form). 

Finally, the maximization over all possible Bk embed- 
ded in Eq. Q is equivalent to find those values {Oopt, 4>opt) 
that disturb the least the spin A when we make a mea- 
sure on B. For Z2-symmetric states we analytically found 
9 opt = = 0, in accordance with Refs. [H1[S]. This 

eventually leads to a simple closed expression for the clas- 
sical correlations: 



C{pr) = SbiniPl) - Sun{P2) ■ 

where Sbmip) is the Shannon entropy and 
l + Q. 1 



Pi 



P2 



(18) 



(19) 



By taking the difference between the mutual informa- 
tion (13) and Eq. (18) for the classical correlations, one 



gets the following simplified expression for the quantum 
discord between any two spins in the XY chain in trans- 
verse field, valid for Z2-symmetric states [5]: 

3 

Q{Pr) = Sm,i{Pi) + Sbin{P2) ^'^^v log2 K ■ (20) 

u=0 



V. QUANTUM DISCORD AT T = 

In this section we analyze quantum correlations both 
in the thermal ground state and in the symmetry broken 
one. In particular we remark differences and similarities 
between them, and highlight the interesting features oc- 
curring at the QCP and at the factorizing field (if not 
specified, the xy anisotropy is set to 7 = 0.7 in every 
picture). 

We start by showing the behavior of quantum discord 
at zero temperature, over a wide range of external field 
values h centered around the critical value /ic, where the 
QPT occurs. In Fig. [T] we plot the quantum discord Qr, 
for both the XY model (main panel, 7 — 0.7) and the 
Ising model (right inset, 7 = 1), in the true ground state 
obtained by means of DMRG simulations (dotted lines 
with symbols) and in the thermal ground state evalu- 
ated analytically (solid lines without symbols). In the 
disordered phase h > he no difference occurs between 




FIG. 1: Quantum discord Qr(h) between two spins at dis- 
tance r in the XY model at 7 = 0.7 (main plot and left inset) 
and at 7 = 1 (right inset), as a function of the field h. Con- 
tinuous lines are for the thermal ground state, while symbols 
denote the symmetry-broken state obtained by adding a small 
symmetry-breaking longitudinal field hx = 10~^ and it was 
computed with DMRG in a chain of L = 400 spins; simula- 
tions were performed by keeping m — 500 states and evaluat- 
ing correlators at the center of the open-bounded chain. For 
7 = 0.7 and at hf ~ 0.714, in the symmetric state afi the 
curves for different values of r intersect, while after breaking 
the symmetry Qr is rigorously zero. At the critical point Qr 
is non analytic, thus signaling the QPT. In the paramagnetic 
phase, there is no symmetry breaking affecting Qr- 



the two state, while in the opposite regime h < two 
different patterns come out. In fact, in the latter case 
the order in the system is really achieved only when the 
symmetry of the system is lost. We achieve this condi- 
tion by introducing a tiny longitudinal field {h^ <^ 1) 
in the DMRG computations that leads to the symme- 
try breaking and gives out the true ground state, where 
quantum correlations are very small as long as h < he. It 
is remarkable that they start to increase once the field is 
tuned immediately upper the factorizing field (where all 
quantum correlations must vanish), eventually reaching 
a cuspid-like maximum at the QCP. On the other hand, 
the quantum discord on the thermal ground state is a 
smooth function with respect to the field. In general it 
depends on the distance r between the two spins, but at 
the factorizing point it gets the same value for any length 
scale, as witnessed by the left inset 

To go deeper in the analysis, let us first focus onhc = 1. 
The QPT is in general marked by a divergent deriva- 
tive of the quantum discord, with respect to the field. 
In particular such divergence is present for any 7 in 
the symmetry-broken state, while on the thermal ground 
state it holds as long as 7 < 1 (see Fig. [2]); for 7 = 1, 
dhQr is finite at he although the d^Qr diverges [5^. This 
divergence suggests that a scaling analysis at the QCP 
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e"^ * \h - Hp 



FIG. 2: Behavior of dhQr in the thermal ground state with 
respect to the field for any type of anisotropy 7. Note that for 
7 = 1 at the QCP there is no more divergence but a cuspid. 
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FIG. 3: Finite-size scaling of dhQ\ for the symmetry-broken 
ground state in the XY model, in proximity of the critical 
point he. The first derivative of the quantum discord is a 
function of L~'^{h — hm) only, and satisfies the scaling ansatz 
dhQ\ ~ L'^ X F[L~'^{h — hm)], where h^. is the renormalized 
critical point at finite size L and uj = 0.472. We found a 
universal behavior he — hm ~ ^-i.28±o.03 .^^j^^j^ respect to 7. 
Inset: raw data of dhQi as a function of the transverse field. 



FIG. 4: Scaling of Qi close to the factorizing field: we found 
an exponential convergence to the thermodynamic limit, with 
a universal behavior according to e~°^(/i — o ~ 1 

[/ij^' denotes the effective factorizing field at size L, while 

5{Qi) = Qi^' — Qi^"*""']- Due to the fast convergence to 
the asymptotic value, already at L ~ 20 differences with the 
thermodynamic limit are comparable with DMRG accuracy. 
Inset: raw data of Qi as a function of h. The cyan line is for 
i = 30 so that, up to numerical precision, the system behaves 
at the thermodynamic limit. 



dence of Qr close to it: 



1 — 7n r 



(21) 



Such behavior is consistent with the expression of corre- 
lation functions close to the factorizing line obtained in 
Ref. [26^ , and here appears to incorporate the effect aris- 
ing from the non vanishing spontaneous magnetization. 
Most remarkably, we found a rather peculiar dependence 
of Qr on the system size, converging to the asymptotic 
value (5^^"*^°°' with an exponential scaling behavior (see 
Fig. 11. 



VI. QUANTUM DISCORD AT FINITE 
TEMPERATURE 



is feasible. In particular in Fig. |3]we show the finite size 
scaling dhQr=i for the symmetry-broken ground state in 
proximity of he- We found that z = v = 1, thus meaning 
that the transition is in the Ising universality class. 

Turning into the factorizing field hf we underline that, 
for the thermal ground state, it is the only value where 
the curves with different r, intersect with each other (see 
up-left inset in Fig. [T]) [21]. Besides this, in the broken 
symmetry state, not only we found that all curves vanish 
m hf, but we numerically estimated the following depen- 



Even if both the QCP and the factorizing field occur at 
T — 0, they influence the physics of the system even if the 
temperature is switched on. Close to he, the physics is 
dictated by the interplay between thermal and quantum 
fluctuations of the order parameter. As we stated before, 
the cross-over temperature T^ross = \h — hc\^ flxes the 
energy scale [TT]. For T <^ Tcross the system is described 
by a quasi-classical theory, while inside the "quantum 
critical region" (T ^ Tcross), it is impossible to distin- 
guish between quantum and thermal effects. Here the 
critical properties arising from the QCP at T = are 



FIG. 5: Quantum discord in the thermal state of the Ising 
model with 7 = 1, as a smooth function of the temperature 
T and of the external field h. 

highly dominating the dynamics of the system; as a mat- 
ter of fact, we expect that quantum correlations display 
some particular pattern as well as they do at he- In fact 
close to hf and at small T, the bipartite entanglement 
remains vanishing in a finite non linear cone in the h — T 
plane |4l[20]. Thermal states, though, are not separable, 
and entanglement is present in a multipartite form |27j . 
In this regime the bipartite entanglement results to be 
non monotonous, and a reentrant swap between parallel 
and antiparallel entanglement is observed [20]. At finite 
temperature, the Z2 symmetry is preserved all over the 
values of h (there is no longer a symmetry broken phase) . 
This means that if the system lies in the ground state at 
T = (symbols lines in Fig. [ij, once the temperature 
is switched on we get a jump of Qr all along the phase 
h < he- After that it behaves as a smooth function de- 
caying with increasing the temperature (Fig. [s]). Such 
discontinuity is also observed in the entanglement, even 
if in that case it is much less pronounced and it occurs 
only for h < hf pT] . 

Let us now analyze how criticality and factorization 
modify the fabric of pure quantum correlations in the h — 
T plane. We start by focusing on the finite-temperature 
scaling of the quantum discord close to the critical point 
he- In the first place we verified the logarithmic scaling 
dhQr\h^ ~ a::ln(r) -|- k along the critical line, h = 1, 
in the h — T plane (see Fig. [9]), where the value of x 
depends on the degree of anisotropy 7. Once x is given 
(for example we found x — 0.065, for r = 1 and 7 = 0.7 
- Fig. [6]), by properly tuning the ratio T/Tcross, where 
Tcross = \h ~ hc\, we verified the scaling ansatz 

duQr^T'^Fi-^Y (22) 

\ ^ cross / 

In particular in Fig. [6] we show how different curves, re- 
lated to different values of T /Tcross, collapse when ap- 
proaching the critical point. Remarkably in the Ising 
case (inset) the scaling is verified as well, even if the 
derivative dhQi is finite at he- To explore the behavior 



T/T 

cross 

FIG. 6: Finite-temperature scaling of the quantum discord for 
the thermal state close to the critical point. The logarithmic 
scaling is verified : along the critical line we found dhQi\ha ~ 
X \n{T)+k, with x — 0.065 for 7 = 0.7. The scaling function F 
shows a data collapse close to the critical point. Inset: same 
analysis for the Ising case (7 = 1); we found an analogous 
scaling behavior with x — —0.0059. 




FIG. 7: Schematic representation of the directional derivative 
on the phase diagram in the h — T plane. It allows to study 
how quantities vary along straight lines coming out from the 
critical point with a slope u = (cos a, sin a). 



of correlations along the h — T plane, we studied how the 
quantum discord varies on the phase diagram just above 
the QCP. In the first place we analyze how the derivative 
with respect to field behaves along the directions fanning 
out from he- In Fig.[7]we sketch a cartoon to describe the 
directional derivative D^Qr — IdxQr sin a -I- daQr cos a] 
we used to describe how Qi varies close to the QCP. 
From the pattern of D^Qi at low temperature (Fig. [S]) 
we see how the presence of the QPT characterizes the 
whole phase diagram. The black vertical line starting 
from the QCP highlights the fact that the quantum dis- 
cord remains constant along the critical line h — 1: in 
a sense, close to such region h ~ 1, quantum correla- 
tions are particularly "rigid" . This explains their robust- 
ness up to finite temperatures, particularly along slopes 
within the quantum critical region. On the other hand, 
out of the quantum critical region, the variation of Qi 
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FIG. 8: Density plot of D^xQi close to he, in the h — T plane. 
The vertical line fanning out from the QCP shows that Qi 
tends to be constant inside the quantum critical region. 
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FIG. 10: Average quantum discord displacement AQr for 
m — 5 fanning out from the factorizing point hf ~ 0.714, 
where all correlations coincide at any length scale r, as evi- 
denced in the left inset of Fig. [T] 
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FIG. 9: Density plot in the h — T plane of drlQi/Ci] close 
to he; along the critical line the ratio Qi/Ci is constant with 
respect to the temperature. The solid straight line ( T = 
\h~ hc\ ) marks the boundary of the quantum critical region. 



is drastically increased. We also point out the pecu- 
liar asymmetric behavior between the two semiclassical 
regions (in the ordered phase D^Qi is generally higher 
than in the paramagnetic phase). Furthermore, to make 
a more accurate analysis we look at the interplay be- 
tween quantum discord and classical correlations. In 
particular we analyze how the ratio Qi/Ci varies with 
the temperature, exploiting the respective sensitivity to 
thermal fluctuations arising at finite temperature. Ac- 
cordingly with the phase diagram related to the QPT, 
a V shaped pattern comes out (see Fig. [9]). In partic- 
ular along the critical line, inside the quantum critical 
region, we found (9t[Qi/Ci] = 0. Then apparently the 
ratio between correlations is constant even though the 
temperature is switched on, as long as the field is tuned 
at the critical value he. Besides this, the whole crossover 
region from a phase to another our is marked as the high- 
est variation in the nature of correlations in the system. 



Finally we concentrate on how the factorizing field af- 
fects the finite-temperature physics of the system. As we 
stressed before, the Z2 symmetry is preserved on the ther- 
mal state. In particular, for the thermal ground state, the 
factorizing field is the unique value where the quantum 
discord is the same at any length scale r. Here we show 
that this feature is present even after the temperature is 
switched on. In Fig. [lO] we propose the quantity 



m(m — 1) ^-^ 



(23) 



as a measure of the robustness of this characteristic at 
non zero temperature. We consider different distances 
between the couple of spins A and B, and take the av- 
erage of the difference between the respective quantum 
discord. Our results strongly support that, for a finite 
range of temperature, this difference is still zero (i.e. the 
quantum discord for different r preserves its scale invari- 
ance only at ft. ~ hf). 

We emphasize that the robustness and sensitivity of 
the quantum discord to non zero temperature, encourage 
the implementation of suitable experiments that could 
give good feedback of our analysis. 



VII. DISCUSSION 

We studied pure quantum correlations quantified by 
the quantum discord Qr in the quantum phases involved 
in a symmetry-breaking QPT. 

In the ordered phase, although Qr results relatively 
small in the symmetry-broken state as compared to the 
thermal ground state, it underlies key features in driving 
both the order-disorder transition across the QPT at he, 
and the correlation transition across the factorizing field 
hf. The critical point is characterized by a non analyt- 
icity of Qr found in the Ising universality class. Close to 
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hf, Qr displays uniquely non trivial properties: in the 
thermal ground state quantum correlations are identical 
at all scales; for the symmetry-broken state the factor- 
ization can be interpreted as a collective reshuffling of 
quantum correlations. We point out that hf marks the 
transition between two "phases" characterized by a dif- 
ferent pattern of entanglement 1151 120] . Accordingly our 
data provide evidence that such a correlation transition 
phenomenon is of collective nature, governed by an expo- 
nential scaling law. We observe that the scaling close to 
hf cannot be algebraic, because the correlation functions 
drop exponentially in the gapped phases. Nonetheless 
we ascribe this specific scaling behavior to the peculiar 
phenomenology of the factorizing phenomenon, which 
goes beyond the generic exponential decay of correla- 
tions that is observed in gapped phases. For finite L 
different ground states for the two parity sectors inter- 
sect [28] . The ground state energy density is diverging for 
all L (such divergence, though, vanishes in the thermo- 
dynamic limit). Indeed we found that the factorization 
occurs without any violation of adiabatic continuity. Ac- 
cordingly, the ground state fidelity J-{h), which can de- 
tect both symmetry breaking and non-symmetry break- 
ing QPT, is a smooth function at /i/ [2^. We remark 
that this can occur without closing a gap and changing 
the symmetry of the system, as a signature of the fact 
that quantum phases and entanglement are more subtle 
than what the symmetry-breaking paradigm says. Such 
a behavior is particularly relevant in the context of QPTs 
involving topologically ordered phases where a QPT con- 
sists in the change of the global pattern of entanglement, 
instead of symmetry [50] . 

We also analyzed the phase diagram at low T. A dis- 



continuity of Qr with T is evidenced in the whole ordered 
phase h < he- We proved that Qr displays universal fea- 
tures, identifying the quantum critical region as the one 
where the quantum discord (relatively to classical corre- 
lations) is frozen out to the T — Q value. In particular in 
each phase the ratio between correlations is stable with 
respect to the temperature, while the highest variations 
develop in the crossover region on the right of he, where 
the system is running out of the critical region into the 
quasi-classical one just above the disordered phase of the 
paramagnet. This type of quantum correlations therefore 
allows a fine structure of the phase diagram, according 
to the behavior of the gap A ^ in the low tempera- 
ture limit T ^ |A|: as a matter of fact, two different 
mechanisms lead to the two corresponding semiclassical 
regimes driven by quantum (A > 0) or thermal (A < 0) 
fluctuations [llj . Finally, a non trivial pattern of quan- 
tum correlations fans out from the factorization of the 
ground state, as well: we identified a finite portion of 
the low-temperature phase diagram where the quantum 
discord is nearly constant at any range, despite of the 
thermal fluctuations. 
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